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Abstract 

The delay difference model was implemented to fit 21 years of brown tiger prawn (Penaeus es¬ 
culentus) catch in Moreton Bay by maximum likelihood to assess the status of this stock. Monte 
Carlo simulations testing of the stock assessment software coded in C++ showed that the model 
could estimate simultaneously natural mortality in addition to catchability, recruitment and 
initial biomasses. Applied to logbooks data collected from 1990 to 2010, this implementation of 
the delay difference provided for the first time an estimate of natural mortality for brown tiger 
prawn in Moreton Bay, equal to 0.031 ± 0.002 week” 1 . This estimate is approximately 30% 
lower than the value of natural mortality (0.045 week” 1 ) used in previous stock assessments of 
this species. 
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1 INTRODUCTION 


1. Introduction 


In fisheries science, estimates of mortality rates are essential to evaluate the status of 
exploited marine resources. The magnitude of natural mortality affects our predictions of 


stock proc 


(Lee et al 


a ctivit y, making this parameter one of the most important to fisheries management 


data (jWand. 


2 01111. Natura l mortality is notoriously difficult, to estimate from catch and effort 


1999 


Zhou et ah 


201ll) . Contemporary fisheries science approaches advocate 


estimation of natural mortality within stock assessment models using catch data to separate 
fishing mortal ity from estimates of total mortality instead of using inaccurate indirect methods 


(Maunder and Pincrj, 


2 014 ). In practice, many researchers are cautious about estimating this 


population dynamic parameter using stock assessments be cause model mis-specifications and 


poor data quality have often led to un-realistic estimates ([Francis . 


201211 . As a consequence, 


the standard practice is to use an externally estimated value of natural mortality as a fixed 
input parameter to fit stock assessment models. 


The stock assessment of brown tiger prawn (Penaeus esculentus ) in Moreton Bay (IKi enzle et al. 


2 0141 1 applied a stan dard approac h by 


ern Prawn Fishery 


estimated by 


Dichmont et al. 


borr owing natural mortality estimates from the North- 


20031). The value of natural mortality (0.045 week 


-i' 


Wan g (119991) using a depletion method has been widely used for over four decades 


without validation. Extending this approach to combi ne est ima tes o f catchability and natural 


mortality for the Northern Prawn Fishery (NPF), ( Zhou et ah 


2 0141 ) reached a substantially 


lower value for natural mortality (approx. 0.03 week x ) suggesting natural mortality might not 
have been as high as previously thought and providing some incentive for further research into 
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1 INTRODUCTION 


this topic. The logic of Wang (119991) ’s approach was to subset catch and effort data using the 
period when recruitment could be assumed to be negligible. But it was based on a single year 
of data from a single stock region and therefore lacked generality. This method would have 
required substantial modifications before being applicable to Moreton Bay fishery to account 
for targeted effort and temperature-dri ven variati ons of availability, two processes identified as 


having major influences on this fishery flKienzle et al 


2 014) . The cost and effort of customizing 


Wang’s method to estimate natural mortality for tiger prawn in Moreton and possibly inte¬ 
grate it to the stock assessment to account for negative correlations between fishing and natural 
mortality might have been greater, if at all possible, than developing the stock assessment itself. 


Therefore, in this paper, we propose to use a different approach to the problem of estimating 


natural mor tality by maximum li 


ence model (Schnu 


since 


Kienzle et al. 


e, 


1985 


celihood for tiger prawn in Moreton Bay. The delay differ- 


Deriso 


1980) used for stock assessment was substantially enhanced 


(20141) to account for additional sources of known variability in this fishery, 


including environmental factors, that improved the predictions of observed variations of catch 
from 1990 to 2010. The delay difference model, implemented in C++, was submitted to an ex¬ 


tensive set of tests to check t 


to expectations (IStroustrupl . 


re qu ality of this program and ensure that it performed according 


20091 ). These tests involving simulations of datasets over parame¬ 


ters ranges consistent with the biology of this species revealed that the model could provide an 
accurate estimate of simulated natural mortality as well as all the other parameters required 
for stock assessment, excluding those pertaining to somatic growth. Therefore we estimated 
natural mortality using real data which provided a realistic value for tiger prawn and suggested 


to use this approach in future assessments. 
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2 MATERIALS AND METHODS 


2. Materials and methods 


2.1. Stock assessment model 


2.1.1. Population dynamics 


A Schnute-Deriso delay-difference model ( Schnutel . 


1985 


Deriso, 


1980 


Hilb or n and W alt ers. 


1992) was used to estimate weekly variations in biomass (B t ) of brown tiger prawns in Moreton 


Bay between 1990 and 2010 


B t — St_i B t _ i + p s t _\ B t _ i — p s t - 1 s t - 2 B t _ 2 — s t - 1 p Wf.~i R t ~ i + uq. R, t . 3 < t < 1092 (1) 


Sex-co mbin ed gro wth parameters (p, Wfc_i and Wk), derived from von Bertalanffy estimates 


(jGri bble and Dredg e, 


19941 ). were fixed in the model (Tab. [[]). This model assumed all prawns 


were fully recruited to the fishery (knife-edged selectivity) at an age of 22 weeks (k = 22), 
weighing 19.5 grams. Wk- 1 , the pre-recruitment weight was interpre ted as a param eter rather 
than the actual weight at age k — 1 and was estimated according to Schnute ( 19851 ). Survival 


(st) varied as a function of a fixed natural mortality rate (M, Tab. CO) and fishing mortality (Ft) 
proportional to effort (F t = qE t ) 


s t = exp [-(M + qEt)} 


( 2 ) 


where catchability (q) was estimated. 
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This model estimated the magnitude of recruitment in each week (R t ) using 1 param¬ 
eter to describe total recruitment in each 21 years (R y ) between 1990 and 2010 and 2-22 
parameters (constant an d year specific /i, k) from the von Mises probability density function 


(iMard ia, and Jupp, 


19991) to allocate a proportion of the total recruitment within each year to 


each 52 weeks 




exp[«; cos(x — n)] 
27 tI 0 (k) 


where Iq{x) is the modified Bessel function of order 0. 


(3) 


2.2. Monte Carlo simulations 


The program that simulated fisheries data was originally develo ped f or quality control pur¬ 


poses, to ensure that the stock assessment software coded in C++ (Ki enzle et al 


be relied on to do what it was supposed to do with an acceptable error rate (iStroustrup . 


2014 


con 


20091). 


ft was designed with Moreton Bay brown tiger prawn fishery in mind, drawing parameters 
values from ranges consistent with this stock (Tab. El) • The dynamics of the fish ery w as sim¬ 


ulated with an age-structured model which alleviated the criticism raised by 


Francis (2012J) 


about validity of simulation testing using the same model to generate datasets and estimate 
parameters. Age (a) varied between 22 and 156 weeks (1 < a < p, p — 135). Each simulated 
dataset consisted of 4 years of data at weekly timesteps ( 1 < t < n, n = 4x 52). 


Individual’s weight at age a > k (where k — 22 weeks is the age at recruitment to the 
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fishery) was calculated according to ISchnutd (119851) Eq. 1.14: 


W t , a = Wt-l,k-l + {Wt,k ~ W t -i,k-l) 


1 ~ P 


1—a—k 


i ~ p 


(4) 


where w t - i,*— i, Wt,k and p where parameters hxed to 17.756, 19.5 and 0.96328 to represent 
brown tiger prawn somatic growth. Weight at age were kept constant between time steps and 
represented using a matrix (W) with dimension n x p having identical rows. 

The value of several parameters were generated at random including catchability, natural 
mortality, von Mises parameters for the recruitment distribution, magnitude of recruitment 
in each year (Tab. [2]). The matrix of effort (E with dimension n x p) had identical columns 
made of a vector of effort (E t ) randomly generated from observed weekly effort in Moreton Bay 
fishery. F is the fishing mortality matrix defined as F — qE where q is a randomly generated 
scalar representing the catchability (Tab. [2]). M is a constant matrix with dimension n x p. 
Natural mortality was generated using a uniform distribution varying between 0.025 and 0.065 
week -1 . 

Recruitment to the population (W,i) was generated using the product of a random yearly 
magnitude of recruitment (Tab. [2]) and distributing across 52 weeks using a von Mises distri¬ 
bution with randomly generated mean and standard deviation (Tab. [2]). 

The matrix of number at age in the population (N) was filled simulating the dynamic of 
n + p — 1 cohorts, applying the declining in abundance along each cohort: 


N i+ ij+i = N itj x exp(—(Mjj + qEij)), fori < i < n - 1 and 1 < j < p - 1 


(5) 
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Hence the matrix of population biomass (B) was calculated using: 


B = N *W 


( 6 ) 


And the vector of biomass for each timcstep: 


b, = j2 b 


(7) 


The matrix of catc h at age (C , with dimension nx p) was calculated using Baranov’s catch 


equation flQuinn and Derisol. 


19991 ) 


C = 


F 


F + M 


B(1 - exp(—(F + M)) 


( 8 ) 


Finally the time series of catch were calculating summing across columns the following 
expression 


c, = '£ c 


(9) 


The time series of catch and effort were altered by ’’white noise” multiplying each element 
by a random factor uniformly distributed between 0.9 and 1.1 to introduce ± 10% error into 
the dataset passed to the stock assessment. 


This simulator was written in R (JR Core Teaml . 


20131) . The code is available from M. 


Kienzle. 
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2.3. Tiger prawns in Moreton Bay 


2.3.1. Water temperature data 


Moreton Bay is a s 


of 36 m (IDavie et al. 


lallow body of sea water, with an average depth of 6.8 m and a maximum 


2 0111) . Since 2002, water temperature has been measured monthly at 85 


stations spread across the Bay by the Ecosystem Health Monitoring Program for the Healthy 
Waterways Program of South East Queensland (http://healthywaterways.org). Station E00527 
is located close to the center of the brown tiger prawn trawling grounds (Fig. [[]). Its bottom 
temperature measurements were used to quantify the emergence behaviour of this benthic 
specie. Maximum daily air tem perat ure were recorded continuously by the Australian Bureau 


of Meteorology (BOM, 


Anonymous) at Cape Moreton Lighthouse, station 040043. A linear 


regression between these two time series of measurements was used to estimate daily bottom 
temperature on the fishing ground. 


2.3.2. Duration of emergence 


Hill (1198511 measured experimentally the effect of temperature (T) on the duration of emer¬ 


gence ( y ) of brown tiger prawn. This author assumed a linear relationship between these 


varia bles whereas he used a sigmoid function for this relationship in crab Scylla serrata (Hill, 


1980 ). Biologically, an asymptotic relationship makes more sense as nocturnal emergence time 


canno t increase infinitely with time. As a consequence, a sigmoid function was fitted to 


Hill 


(ll985l Vs experimental data 


y(T) 


a 

1 + exp(& - cT) 


( 10 ) 


in place of a linear relationship using the non-linear regression functions in R (1R, Core Team! . 
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20131 ). Goodness of fit was assessed with a y 2 statistics (IBrandtUlOOSI) . An index of availability 
(y(T) varying between 0 and 1 at a given temperature, in °C, was calculated using the predicted 
duration of emergence ( y(T )) divided by the asymptote of the sigmoid model 


7 (T) = 


y(T) 


(ii) 


This in dex of availability was used as a multiplier to effort in the stock assessment model 


(JKienzle et ah 


2 0141) to test the assumption that emergence time affected catchability. 


2.3.3. Refining the model to estimate natural mortality 

The present work extended the delay difference model presented in 


Kienzle et ah 


(2014). 


For consistency, model numbers were kept identical and the reader is referred to the previous 
study for details on model 1 to 7. 


The first area investigated concerned the prawn emergence model: model 8 replaced the 
linear model in model 4 by a sigmoid relationship. Model 9 used estimated daily water temper¬ 
ature rather than averaged monthly temperature off Cape Moreton. Model 10 was identical to 
model 9 but estimated natural mortality rather than fixing it to 0.045 week” 1 . Finally model 
11 allowed for timing of recruitment to vary from year to year. 
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3. Results 

3.1. Monte Carlo simulations 

Monte Carlo simulations showed that the implementation of the delay difference in C++ was 
capable of estimating natural mortality, catchability, von Mises parameters mean and standard 
deviation and initial biomasses (Fig. [2]). Simulated values for natural mortality over the range 
0.025-0.065 week -1 (Tab. [ 2 ]) containing realistic values for tiger prawn could be estimated with 
less than 31% discrepancy from simulated values in all cases; less than 15% discrepancy in 95% 
of cases and less than 4% discrepancy in 50% of cases. By comparison, catchability estimates 
were up to 67% different than their simulated value with 95% of cases within 18% discrepancy 
and 50% within 2.5%. 


The AIC was used to measure the benefit of estimating natural mortality: the larger the 
discrepancy between fixed natural mortality (0.045) and simulated values, the larger the im¬ 
provement in AIC obtained by fitting the model that estimated natural mortality compared to 


the one that did not (Fig 3(a)). 


Fixing natural mortality in the stock assessment to a specific value influenced estimates of 
catchability: fixed over-estimates of natural mortality produced under-estimates of catchability 


and vice-versa (Fig 3(b)). The negative correlation between these 2 variables was p = —0. 


3.2. Daily temperature estimates and trends 

The maximum daily air temperature, collected by the BOM at Cape Moreton, decreased 
slightly from 1957 to the mid 1980s and increased in average by 1°C from the mid 1980s to the 
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end of 2011 (Fig. 4(a)). A linear regression estimated air temperature increased in average by 


0.180 ± 0.014 °C per decade over the entire period. 


An analysis of variance (ANOVA, Tab. |3|) showed that variables most influencing sea tem¬ 
perature in Moreton Bay were, by decreasing order of importance: time; location and depth. 
The vertical gradient in temperature was very small because layers of water mixed well in 
this shallow Bay. The difference in water temperature between surface and deeper layer was 
0.0162 ± 0.012 °C.m _1 , approx. 0.1°C at the average depth of 7 m. 


A linear regression between 85 matching pairs of sea and air temperature measurements (as 
a dependent variable) was statistically significant (F-statistic: 354.8 on 1 and 83 DF, p-value: 


< 2.2 10 16 , adjusted R 2 = 0.81, Fig. 4(b)). The intercept did not differ significantly from zero 
and the estimate of slope was equal to 0.89112 ± 0.04731 (Tab. SI). The range of estimated 
daily bottom temperature inside the Bay, at station E00527, between 1988 and 2011 was larger 
(11.2-32.2°C) than ocean water temperature measured off Cape Moreton (20.2-26.1°C). Es¬ 


timated daily sea temperature in the Bay have increased between 1990 and 2010 (Fig. 4(c)) 
by an estimated 0.6°C using yearly averages in each of the two decades 1990-1999 and 2000-2011. 


3.3. Model for emergence duration 

The sigmoid model fitted well the experimental data on emergence time (% 2 9 df = 6.81, P 


= 0.34, a = 444.5 ±61.11, b = 9.650 ±2.394 and c = 0.4702 ± 0.1306, Fig 4(d)). The estimated 


index of availability (y(T)) varied periodically over 12 months, reached its lowest values in 
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southern hemisphere winter and highest in summer (Fig. 4(e)). Its amplitude encompassed 


differences in availability of 90% between the warmest summer and the coldest winter. 

3-4- Stock assessment model 

Allowing the timing of recruitment to vary from year to year (model 11) provided the largest 
stock assessment model improvement (AAIC = -272.21, Tab. [2D compared to the second best 


model (model 10). Mean recruitment was estimated to vary from Jan to March (Fig. 3(d)) 


Previous work had found that temperature-driven behavioural response that changes catch- 
ability seasonally provided large improvements in fitting the delay difference to catch data 
(models 5-6 and models 3-4, Tab. 2]). The sigmoid model for emergence duration was slightly 
better than the linear model (compare models 8-3 and 4-3, AAIC = -168.97 and AAIC = - 
164.71 resp.). On the other hand, using daily estimates of bottom temperature provided a slight 
improvement compare to using average monthly temperatures (compare model 9-4 and 9-8, 
AAIC = -49.86 and AAIC = -45.6) but the introduction of larger variations in the availability 
index decreased catchability estimates by nearly half, for both targeted and non-targeted effort. 


Estimating natural mortality improved further the fit of the delay difference (model 10 in 
Tab. E]). In this case, natural mortality was estimated to be equal to 0.035 ± 0.002, about 25% 
lower than its assumed value of 0.045 wk” 1 in other model. This additional degrees of freedom 


in the delay difference also narrowed the distribution of the recruitment (Fig. 3(d)), probably 


because higher survival meant reduced recruitment magnitude required to explain observed 
catches. Finally, the model that best fitted the data (model 11) estimated natural mortality 
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equal to 0.031 ± 0.002. 


4. Discussion 


Monte Carlo simulations showed that natural mortality could be estimated within this stock 
assessment model. The model that best fitted brown tiger prawn data from Moreton Bay es¬ 
timated natural mortality to be equal to 0.031 ± 0.002 week -1 . This estima te is consistent 


with that for common banana prawns (Penaeus merguiensis ) obtained by 


Zhou et ah 


(2014) 


using a depletion method. Both these estimates are approximately 30% lower than the value of 


natural mortality common 


on the analysis by 


y used in brown tiger prawn stock assessment (0.045 week x ) based 


Wan g ( 19991 ) who concluc 


and no more than 0.065 week 1 . 


Zhou et ah 


ed that natural mortality is at least 0.03 week 


-i 


(20 14) attributed this discrepancy to differences 


in behaviour between tiger and banana prawns: the latter forming large aggregations which 
could improve their chance of surviving predators attacks. An alternative explanation of the 
difference between natural mortality estimates might be that WanJ ( 1999 ) used dat a cove ring a 


smalle r area and shorter period of time, October to November in a single year, than 


Zhou et al. 


(120141 ). evidencing spatial or temporal variations of natural mortality rates of prawns across 


their habitat due for example to increased predation in spring. Finally, this discrepancy could 


also have arisen from the difficult y o f 


multaneously reported by 


estimating both catchability and natural mortality si- 


Wan g (1l999h . who resorted to fix catchability: in this situation, an 


under-estimate of catchability could have led to over-estimating natural mortality. At present 
the information about natural mortality estimates is too sketchy to draw firm conclusions. 


Further research focused on developing additional independent estimates of natural mortality 
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should help identify the cause(s) of such discrepancy. 


The conclusion that Moreton Bay brown tiger prawn dataset seems to allow estimation of 
natural mortality within its stock assessment is difficult to generalize because a big impediment 
to progressing fisheries science is the availability of good quality data. That Moreton Bay offers, 
as far as we can see, a reliable logbooks dataset for tiger prawn that is easily fitted by the delay 
difference model is a meagre consolation against the backdrop of many fisheries which data qual¬ 
ity varies from questionable to unusable. Stock assessments scientists often resort to external 
parameter estimates when data can’t provide a reasonable estimate of important aspects of the 
dynamic of the stock. To a large extent, all assessment models are mis-specihed when dealing 
with real data: they are simplification of a complex reality described only partially by available 
data. That is, they include many assumptions that we know are false in the real world, but that 
are needed to make estimation tractable ( e.g ., it is often assumed that neither fishery selectivity 
nor natural mortality vary from year to year). We cannot avoid mis-specihcation because the 


real wor’ 


models (iFrancisl . 


d is co mplex and heterogeneous in ways that we can never expect to reproduce in our 


201 21). ’’Essentially, all models are wrong but some are useful” (George E.P. 


Box). There are no doubts that the present model is also mis-specihed: for example, trawl 
selectivity have surely changed throughout the period investigated in an unknown (and impos¬ 
sible to know) fashion, yet this implementation of the delay difference model assumes constant 
knife-edge selectivity fixed throughout the time-series. Ignoring this change in the fishery has 
certainly had an effect on natural mortality estimates of an unknown magnitude and direction. 
But the fundamental question raised by this paper is whether to use mortality estimates that 
ignore all we have learned about the dynamic of this fishery and use, without questioning, a 
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value borrowed from another fishery. Or to allow the stock assessment model to estimate it 
simultaneously and work at removing as many mis-specihcations as we can in order to converge 
toward a model that provides the best values to parameterize the dynamic of this fishery. 


Natural mortality is regarded as one of the most influential quantities in fisheries stock 
assessment and management. The magnitude of natural mortality relates directly to the pro¬ 


ductivity of the st ock, t 


points (Lee et al 


2011 


le yields that can be obtained, optimal exploitation rates, and reference 
Since information on total mortality is contained in the catch data, 


any error in the chosen value of natural mortality, within limits, is compensated by an error 
in opposite direction on fishing mortality rates: an over-estimate of natural mortality leads to 
under-estimate catchability and the resulting F will be under-estimated and stock biomasses 
over-estimated. The stock may appear to be better than it truly is, and fishing mortality 
appears less of a concern than it should be. For some data poor fisheries, the target fishing 
mortality is often linked to M (e.g., setting F msy = M), in such cases over-estimating M will 
cause stock being over-fished. Reference points estimates, such as maximum sustainable yield, 
were shown to be robust to natural mortality mis-specihcations in a simple age-structure stock 


assessment flCl arki . 


19991) . a modest over-estimate of the natural mortality can result in quota 


recommendations that represent very high real harvest rates. This author recommended to use 
a conservative (low) estimate of natural mortality for management purposes because underes¬ 
timating natural mortality rate carry little penalty in terms of long term yield. Investigating 
the effects of using erroneous natural mortality rate in the delay-difference model on reference 


points was beyond the scope of this work. But such study would clarify whether using these 


most recent, smaller, estimates of natural mortality has a significant impact on harvest deci- 
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sions taken to manage prawn stocks. 
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Figure 1: Location of the sampling site in Moreton Bay. 
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Figure 2: Results from 100 simulations testing the capability of the delay difference model coded in C++ to 
estimate parameters. Each panel compares simulated values (x-axis) to estimated values (y-axis). The vertical 
bars represent 1 standard deviation. Panel (a) natural mortality; (b) catchability; (c) von Mises distribution 
mean; (d) von Mises distribution standard deviation and (e) initial biomasses {B\ and i? 2 ). The dash line 
represents y = x. 
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Figure 3: Effects of fixing natural mortality in the delay difference model: (a) differences in AIC between model 
fixing natural mortality and model estimating natural mortality (y-axis) as a function of the simulated value of 
natural mortality (larger values indicated better fit by estimating natural mortality); (b) effect of fixing natural 
mortality on catchability estimates and (c) distribution of natural mortality estimates for each value fixed in 
the simulation and (d) Von Mises distributions using parametrization estimated with different delay difference 
models. 
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Figure 4: Variations of temperature and its effects on availability: (a) time series of maximum air temperature 
collected by the Australian Bureau of Meteorology at Cape Moreton. The smooth curve represents the trend 
in the data; (b) linear regression between sea temperature collected at the center of the Bay (station E00527 
on y-axis) and maximum daily air temperature recorded at Cape Moreton (x-axis); (c) estimated temperature 
at location E00527; (d) expe rimental data on the effect of temperature on the duration of emergence of brown 
tiger prawn from fUilil ( 19851) . The grey vertical bars represent an estimated variability of measurments. The 
black contin uous li ne re presents the fit of a sigmoid function. The grey dotted li nes represen ts the linear model 
presented bv iHilll ( 1985h and used in previous stock assessment models ( Kienzle et all [2014 ) and (e) estimated 
weekly index of availability. 
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Tables 


Parameter 

Value 

Reference 

P 

0.963 

Gribble and Dredge 119941 

Wk-l 

17.8 grams 

based on the method of Schnute 119851 

Wk 

19.5 grams 

based on the method of Schnute 119851 

M 

0.045 week" 1 

Dichmont et al. (2003) 


Table 1: Values of parameters fixed in the delay difference model. 


Variable type 

Distribution 

Parameters 

recruitment 

uniform 

min=le6, max=le8 

von Mises mean 

uniform 

min=—7T + 0.2, max = n — 0.2 

von Mises sd 

uniform 

min=l, max=50 

natural mortality 

uniform 

min=0.025, max=0.065 

catchability 

uniform 

min=le-4, max=le-3 

fishing effort 

fixed 



Table 2: Distribution and range of value taken by different variables in simulations. 



Df 

Sum Sq 

Mean Sq 

F value 

Pr(>F) 

Location 

84 

4366.38 

51.98 

42.90 

0.0000 

Depth 

1 

445.40 

445.40 

367.61 

0.0000 

Year:Month 

367 

606282.41 

1652.00 

1363.47 

0.0000 

Residuals 

58585 

70982.03 

1.21 




Table 3: Analyis of variance of sea water temperature measurements. 
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Estimate 

Std. Error 

t value 

Pr(> t ) 

(Intercept) 

0.0486 

1.0907 

0.04 

0.9645 

X 

0.9081 

0.0448 

20.26 

0.0000 


Table 4: Linear regression between sea water temperature in Moreton Bay and maximum air temperature at 
Cape Moreton. 






Model 

P 

AIC 

Catchability 



Recruitment distribution 

Initial biomasses 



Nat 

ural mortality 




qi (x 10~ 4 ) 

Q2 (Xl0- 4 ) 

P 

K 

Bi (xlO 5 ) 

B 2 (xlO 5 ) 

a 


M 

11 

49 

6743.71 

3.45 ± 0.532 

1.424 ± 

0.22 

0.668 ± 0.215 

4.955 ± 0.591 

1.012 ± 0.015 

0.978 

± 

0.011 

5.073 ± 0.114 

0 

031 ± 0.002 

10 

29 

7015.92 

2.245 ± 0.377 

0.839 

± 

0.149 

0.525 ± 0.03 

4.734 ± 0.586 

2.173 ± 0.682 

2.038 

± 

0.615 

5.853 ± 0.127 

0 

035 ± 0.002 

9 

28 

7040.14 

2.614 ± 0.801 

1.052 

± 

0.285 

0.61 ± 0.061 

2.819 ± 0.391 

2.078 ± 0.82 

1.923 

± 

0.741 

5.924 ± 0.132 



7 

28 

7071.04 

4.089 ± 0.425 

1.858 

± 

0.255 

0.642 ± 0.036 

2.022 ± 0.088 

0.457 ± 0.011 

0.447 

± 

0.012 

6.008 ± 0.13 



8 

28 

7085.74 

4.761 ± 0.485 

1.898 

± 

0.247 

0.626 ± 0.035 

2.053 ± 0.09 

0.375 ± 0.15 

0.372 

± 

0.131 

6.049 ± 0.133 


a 

4 

28 

7090.00 

4.164 ± 0.44 

1.713 

± 

0.236 

0.649 ± 0.036 

2.047 ± 0.09 

0.451 ± 0.012 

0.444 

± 

0.013 

6.06 ± 0.133 


i-s 

P 

6 

27 

7162.90 

1.612 ± 0.215 




0.513 ± 0.022 

5.29 ± 0.582 

1.025 ± 0.342 

0.979 

± 

0.301 

6.272 ± 0.132 


SP 

3 

28 

7254.71 

4.268 ± 0.357 

1.888 

± 

0.222 

0.499 ± 0.032 

2.358 ± 0.128 

0.148 ± 0.079 

0.163 

± 

0.068 

6.535 ± 0.144 



2 

28 

7265.15 

5.034 ± 0.42 

1.883 

± 

0.228 

0.471 ± 0.032 

2.254 ± 0.118 

0.128 ± 0.068 

0.142 

± 

0.059 

6.566 ± 0.141 



5 

27 

7332.61 

3.243 ± 0.283 




0.513 ± 0.026 

3.386 ± 0.269 

0.182 ± 0.09 

0.189 

± 

0.077 

6.778 ± 0.147 



1 

27 

7369.77 

3.463 ± 0.299 




0.486 ± 0.023 

3.612 ± 0.31 

0.17 ± 0.087 

0.177 

± 

0.075 

6.894 ± 0.154 




Table 5: Comparison between several models fitted to catch using a negative log-likelihood function. Results were ordered by increasing values of 
Akaike Information Criterion (AIC) from top to bottom. 
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